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ABSTRACT 

Recent results on Dalitz analysis of three-pseudoscalar decays are discussed in the 
context of probing charm hadronic-decay mechanisms: the role of FSI effects, which 
create phase shifts between the interfering resonant channels, can be studied in 
the different decay modes and the annihilation contribution measured in the charm 
sector through the Df — > 7r + 7r~7r + decay. 



1 Introduction 



The analysis of substructures should be considered in the context of probing hadronic- 
decay dynamics. The quality of data and the availability of high-statistics samples 
have already shed a great deal of new experimental light on the non-leptonic de- 
cays of charm particles. Indeed, a full lifetime-hierarchy pattern has now been 
experimentally established E' H), posing severe constraints on theoretical models 
and pointing to hadronic-sector corrections. The importance of final-state inter- 
actions (FSI), in which the two decay daughters undergo strong rescattering after 
their initial formation, can be gauged from the isospin-amplitude interference. The 
factorization model predictions can be seriously distorted by FSI; the example of 
D° — > K + K~ /7r+7r~ ratio, expected to be 1.4 (BSW ED) and measured as ~ 2.5 ID, 



Table 1: Isospin Decompositions. 



Mode 


Ratio of amplitudes 


6 = §1 — 5p 


Ktt 




^1/2/^3/2 


= 3.99 ±0.25 


86° ± 8° 


K*tt 




^1/2/^3/2 


= 5.14 ±0.54 


101° ± 14° 


Kp 




^1/2/^3/2 


= 3.51 ±0.75 


0° ± 40° 


K*p 




^1/2/^3/2 


= 5.13 ± 1.97 


42° ± 48° 


KK 


\M/\M -- 


= 0.61 ±0.10 


47° ± 10° 


7171 


IM/IM -- 


= 0.72 ±0.14 


82° ± 9° 



underscores their important role. In table [I] a list of phase shifts and isospin ampli- 
tudes is reported: the phase shifts between different isospin amplitudes are seen 
to be often approximately 90°, indicating the general importance of FSI. 
Recently, amplitude analysis of non-leptonic decays has emerged as an excellent tool 
for studying charm hadron dynamics. In fact, deviations from flat Dalitz plots in- 
dicate non-constant dynamics and a knowledge of the quantum-mechanical decay 
amplitudes allows one to properly account for interference effects when calculat- 
ing branching ratios. In addition, such an analysis provides new probes of FSI by 
measurement of the phase shifts between the interfering amplitudes for the various 
resonant channels. A problem, still unsolved in charm decay, is the reliable esti- 
mate of non-spectator contributions; Df — > 7r + 7r~7r + is the best candidate to occur 
through an annihilation diagram since the annihilation of the two initial quarks 
is Cabibbo favoured and not suppressed as in the D + . Nevertheless, the presence 
of resonant channels with an ss quark content would suggest spectator processes 
rather than W annihilation, making a Dalitz-plot analysis necessary to disentan- 
gle direct three-body decays from resonances. Many experiments have focused on 
the substructure analysis; among the simplest decays to analyze are those of the 
ground-state charm mesons into three pseudoscalars: 

D -> Knn (Mark III D , E691 0) , Argus D , E687 D ) 

D -> KKtx (E691 ® , E687 0) ) 

D -> TXTXTx (E691 0), E687 [preliminary]). 

2 KKtt Dalitz plots 

The Df, D + — K + K~tt + [] Dalitz plots in m 2 KK , m 2 KlT , obtained by E687 and shown 
in fig. 0, are particularly instructive. The Dalitz is very highly dominated by 

1 Through this paper the charge conjugate state is always implied. 



the 07r + and K*°K + channels while the D + has a significant additional contribution 
appearing as a uniform event intensity. Depopulation of the central region of the 
<p and K*° is due to a node in the decay function, reflecting angular-momentum 
conservation; the cos Oheiidty (fig- @) dependence is due to the vector nature of the 
mesons (0, K*°) decaying into two pseudoscalars. 



Figure 1: KKn mass and Dalitz plots of the D + and D s mass region. 




Figure 2: Helicity angle definition. 



The Dalitz plots are fit with a coherent sum of resonances: quasi-two-body channels 
of the type 

D — > r + c 

U a + b, 



which are described by a decay function (eq. [I] and fig. |3|), 

M = F D F r x \c\ J \ a\ J Pj (cos QQ x BW(m ab ). (1) 
The matrix element is the product of 

• two vertex form factors (Blatt-Weisskopf momentum-dependent factors), 

• an order J Legendre polynomial representing the decay angular wave function, 

• a relativistic Breit-Wigner (BW). 




q 2 -(m -iT/2) 2 




Figure 3: D — > KKtc decay diagram. 

The total decay amplitude consists of various functions of the form (|l|), each mul- 
tiplied by a complex factor a r e lSr ; in the simpler case of the Df, where only two 
contributions dominate (07r + and K*°K + ) the total amplitude is 



K + K-^ 



a R * e i5 K*°M (n + K~K + \K(892) 



+ a^M{K + K--R- 



(2) 



where the modulus a measures the importance of the channel in the decay and S 
is the phase. With the factorization hypothesis, the bare amplitudes will be real. 
Phases are then induced by FSI as consequence of Watson's theorem, which states 
that the observed amplitudes are related to the bare amplitudes by the square root 
of a strong-interaction S-matrix, describing the hadron-hadron scattering, 

' r]e 2iSl — r] 2 e i( - Sl+s ^ \ 1 ^ 2 ' 



' obs 



rje 



2,6-, 



a F , \ 

^ hare 



(3) 



through the strong- interaction phase shifts ($1,2) and the elasticity parameter (77). 
While FSI can be measured in two-body decays through the interference of isospin 



amplitudes, in three-body decays they can be inferred from the interference of res- 
onant channels, i.e., from the phase shifts. In practice, the phase of one resonance 
channel is factored out and hence relative phases are measured from experimental 
fits to A* A. Visual evidence for FSI effects is possible through the study of the 
behaviour of the K*° lobes in D + — > K + K~tt + (fig. [I]). The asymmetry, not present 
in the £)+ case, can be interpreted as an interference effect (fig. |j) with a nearly 
constant amplitude (a broad scalar), which may be written as cos 5 + i sin 8. The 
interference of these two amplitudes will produce a contribution to the intensity of 
the form 




1 1.5 2 2.5 3 1 1.5 2 2.5 3 



m 2 K+K _(GeV/c 2 f m 2 K+K _(GeV/c 2 f 

Figure 4: The interference pattern of D + — > KK-k in the K*° band. 

As one moves along the K*° band in the direction of increasing m 2 KK , cos 8 changes 
from +1 to —1 and both terms in eq. [| switch sign. Along a line of constant 
m? KK the first term, which dominates for a real amplitude (5 ~ 0°), switches sign 
on passing through the resonance pole (thus cancelling the interference), while the 
second interference term, which dominates for a relatively imaginary amplitude does 
not change sign under this motion. Thus, the patterns for real and imaginary phases 
are quite different; the interference present in the D + more closely resembles the 
pattern for a relatively imaginary amplitude, 5 ~ 90°. A fully coherent Dalitz 
analysis of D + — > K + K~n + confirms this supposition (see table 2). 



2.1 Speculation on the universality of Dalitz phases 

Since the phase factors in the D + and D+ both originate from FSI, one is tempted to 
speculate that they might be the same. In fact, if the D + and were degenerate 
in mass, the S"-matrix describing the rescattering for 07r + and K*°K would be the 
same: 

/ rye 2 ^ 1 - r]' 2 e i( - Sl+S2) 

~ \ iy/l - n 2 e i{&1+52) r]e 2i52 

For elastic FSI (77 = 1) it is easy to see that the observed amplitudes would then 
pick up a relative phase shift given by 5\ — 62, which is common to and D + 
in the limit of degenerate masses. In this case, defining R = (a^/a^* )p% s . and 

(?)"(? A f(?) 

and thus 5^ — <%* = 01$ — flf,o would be equal for both D + and Df. If r] < 1, 
the observed phases would depend on the ratio of the (f> and K*° amplitudes, which 
are different in the and D + and thus 5$ — 5j^ t0 (D + ) ^ 8$ — Sj^, (D^). In the 
particular case of total inelasticity, rj = 0, 

(*H:;fm^(:;)m « 

and 5^ - 8 K * = 90° - 2 tan" 1 R . 

Thus, the universality of Dalitz amplitude phases might provide a clue to the elas- 
ticity of FSI. Table 2 compares the E687 final results for D + and 0) obtained 
through a fully coherent analysis, the total fit amplitude containing various K + K~ 
and K~tt + resonances plus a flat non-resonant term. A few comments may be made 
on these results. The D + consists of nearly equal contributions of K(892)*°K + , (pir 
and the broad scalar K(1A30)*°K while the is strongly dominated by just K*°K + 
and 4>tt + plus a small contribution from / (980) and /j(1710), which sum to 14%. 
Both the D + and can be fit by entirely quasi-two-body processes without the 
inclusion of a non-resonant contribution. It is interesting to note that both charm 
states have relatively real phases between the dominant K*°K + and 07r + channels, 
sin [5^ — <%„o) ~ 0, which indicate absent or cancelling FSI phase shifts. In con- 
trast, for the D + , sin (<%( 1430 )»o — <%(892)*°) ~ 0.94, pointing to strong FSI effects, as 
anticipated in the discussion of the asymmetric K*° lobes, with the i^(1430)*° play- 
ing the role of the nearly constant scalar resonance. Unfortunately, the large error, 




Table 2: Comparison of D + and Df — > K + K ir 



Parameter 


D+ 


Df 




0° (fixed) 


0° (fixed) K* 




-159 ± 8 ± 11° 


178 ± 20 ± 24° 


<%(1430)*°A" 


70 ± 7 ± 4° 


152 ± 40 ± 39° 


^/o(980)tt 




159 ± 22 ± 16° 






110 ± 20 ± 17° 


/i?(892)*°A" 


0.301 ± 0.020 ± 0.025 


0.478 ± 0.046 ± 0.040 


f(j)TT 


0.292 ± 0.031 ± 0.030 


0.396 ± 0.033 ± 0.047 


fK(1430)*°K 


0.370 ± 0.035 ± 0.018 


0.093 ± 0.032 ± 0.032 


//o(980)tt 




0.110 ± 0.035 ± 0.026 


/ fj(1710)K 




0.034 ± 0.023 ± 0.035 



(t(AS) = 56°, on the 7^(1430)*°^ relative phase for the Df precludes a meaningful 
test of the conjecture of a universal phase shift for this channel. 

3 E687 preliminary results on Df, D + — > 7r + 7r _ 7r + 

The observation of a non-resonant three-pion decay of the Df would represent a clear 
signature in favour of non-spectator processes in the charm sector. To determine 
the rate of Df — > 7r + 7r~7r + due to annihilation, the contributions from the resonant 
decays that most likely originate from the ss quark through spectator-diagrams have 
to be removed. The power of the E687 analysis is the air-gap cut, which imposes 
the requirement that the secondary vertex lies in the air-gap region downstream of 
the Be interaction target and upstream of any spectrometer elements. It signifi- 
cantly reduces the non-charm background and increases the signal-to-noise ratio by 
about a factor of four in both the D + and Df cases. In fig. |5| the invariant mass 
plots are shown with and without the air-gap cut. The large Ktnr contamination, 
with the K misidentified as a 7T, is removed by requiring the candidate mass to be 
incompatible with the D + — > Kirn hypothesis. The Dalitz plots for the D + and 
Df signal regions, defined within ±2cr from the mass peak, are shown in fig. |6|. 
Owing to the presence of two identical particles, the 7r + 7r~ amplitudes have to be 

2 The fractions f r are known as decay fractions and represent the ratio of the integrated Dalitz 
intensity for a single resonance r divided by the intensity with all contributions present 

/ dml^ml^aie^Miabc^Yia^Miab^r)) 
decay raction, - j dm^dmU^ M(abc\r))* M(abc\r)) ' 



Bose symmetrized; the representation chosen is that of a folded Dalitz plot where 
the lower and the higher value of the two possible mi + combinations are plotted 
on the abscissa and ordinate respectively. 




7C7C7C mass spectrum tctctc mass spectrum with air gap cut 

Figure 5: mm mass spectrum a) without and b) with the air-gap cut. 

Following the model already explained, the total amplitude is assumed to consist of 
a flat uniform term representing the non-resonant contribution plus a sum of func- 
tions representing the intermediate strong resonances and decay angular-momentum 
conservation. The fit parameters are then the amplitude coefficients a, and phases 
8i. All the resonances decaying into 7r + 7r~ with sizeable branching fractions have 
been considered in the fit function. Fig. [7] shows the two mass projections of the 
Dalitz plots together with their sum. The resulting fit fractions and the relative 
phases are listed in table |3|. 

The D + appears to be largely dominated by the non-resonant and pvr + , with marginal 
contributions from the scalar /o(980) and the spin-2 resonance /2(1270). A signature 
for the p(770) is the peak in the mass projections around 0.6 (GeV/c 2 ) 2 and for the 
/o(980) the bump around 1.0 (GeV/c 2 ) 2 . The Dalitz plot, in contrast to the D + , 
is more structured. Even a visual inspection reveals the / (980) band, previously 
observed by E691 0), and an accumulation of events at the sharp corner on the 
right side of the Dalitz plot, which could be attributed to the third helicity lobe of 
the / 2 (1270) plus an accumulation near the contour on the left side corresponding to 



Table 3: D + and — > 7r + 7r 7r + comparison. 



Parameter 




Df 




0° (fixed) 


235 ± 22° 


U f)7r 


27 ± 14° 


53 ± 44° 


^/ 2 (1270)7r 


207 ± 17° 


100 ± 18° 


0/ (980)tt 


197 ± 28° 


0° (fixed) 


^/ 5 (1500)7T 




234 ± 15° 


fnr 


0.589 ± 0.105 


0.121 ± 0.115 




0.289 ± 0.055 


0.023 ± 0.027 


// 2 (1270)vr 


0.052 ± 0.034 


0.123 ± 0.056 


//o(980)tt 


0.027 ± 0.031 


1.074 ± 0.140 


// S (1500)7T 




0.274 ± 0.114 



the first helicity-lobe region of the /2(1270). The non-resonant contribution appears 
to be marginal because of the absence of population in the upper corner of the Dalitz 
plot. 

By way of contrast, a uniform-density horizontal band around 2.0 (GeV/c 2 ) 2 is 
manifest, whose structure suggests the presence of a scalar resonance in 7r + 7r~ mass 
around 1.5 (GeV/c 2 ). Given the major uncertainty surrounding dipion scalar reso- 
nances in this mass region, the experiment decided to use its own data to find the 
parameters for a single state, which best reproduces the Dalitz plot. Several fits of 
Df Dalitz plot were performed by varying the mass (1200 MeV - 1600 MeV) and 
the width (25 MeV - 400 MeV) of this state, which was parametrized as a Breit- 
Wigner resonance. Fits were performed taking into account all possible additional 
contributions from a non-resonant component and the following well-established res- 
onances, / (980), / 2 (1270) and p(770). A clear maximum of the likelihood function 
is found at m = 1.475 GeV/c 2 and F = 100 MeV/c 2 . The parameters of this state, 
which is denoted as 5(1500), are remarkably consistent with the /o(1500) entry of 
PDG96.0 Evidence for the non-resonant channel, as well as the pn + contribution, 
would have important implications for the theory; while the latter decay is expected 
to be heavily suppressed in the factorization model assuming isospin symmetry 0) , 
the evidence of the non-resonant channel would indicate the presence of an initial- 
quark annihilation process. However, in contrast with the p(770)7r + case, assessing 
the presence of the non-resonant component in a decay is a very difficult task, which 

3 Although the mass and width of this state are in excellent agreement with the /o(1500) entry 
in PDG96, it is to be emphasized that several interfering resonances could equally well describe 
the D+ Dalitz plot. 



would at least require a high-statistics sample since a coherent sum of wide reso- 
nances could easily mimic an almost flat contribution. A comparison of the D + 
and — ► 7r + 7r~7r + phase shifts (table 0) suggests a different and presumably more 
complex dynamics with respect to the K + K~tt + , with a different role of the FSI for 
D + and D£ in this decay channel. 




Figure 6: Folded Dalitz plots for a) D + and b) D+ — > ixixix . 



4 Kittt Dalitz plots 

The Ktttc final state is the most studied so far: many experiments have focused 
on it and results are published for D + — > fC~7r + 7r + , D° — > fT|?7r + 7r~ and D° — > 
i^ _ 7r + 7r°. A certainly striking conclusion reached by all groups is the presence of a 
large non-resonant component for D + — > K^7i + 7i + , which is unique to this decay, 
the majority of charm-meson decays being dominated by quasi-two-body processes. 
Furthermore, amplitude interference is significant in all three modes; the sum of the 
decay fractions, reported in table [|, greatly exceeds unity EP . 

An analysis of the phase shifts is interesting also in this case: in table [5] the E687 
| sin (5 — <%*„.)| is reported for the three modes. 

It is interesting to note that the sinA<5 values are consistent row by row: inelastic 
FSI or isospin phase shift could break this. One might also be tempted to affirm the 



Table 4: Sum of the Kirir final state decay fractions. 





K 7T + 7T + 






E/r 


1.47 


1.27 


1.18 



Table 5: | sin (5 — o~x* n )\ for the three D + — > fC7r7r modes. 





K 7T + 7T + 






Kp 




0.69 ± 0.1 


0.31 ± 0.2 


K*(U30)7r 


0.26 ± 0.09 


0.24 ± 0.23 




nr 


0.74 ± 0.03 




0.97 ± 0.25 



presence of large FSI for K*tt relative to the non-resonant channel, but the issue here 
is to understand the precise nature of the non-resonant contribution: i.e., whether 
it is a single process or a mixture of broad resonances. In addition, although the 
amplitude model qualitatively reproduces many features of the data, statistically 
significant discrepancies are observed in some of the fits. This may be suggestive 
of strong-interaction dynamics not contained in the model, such as the presence of 
new, undiscovered wide resonances or a non-uniform, non-resonant amplitude with 
a possibly varying phase. 

5 Conclusions 

Dalitz-plot analysis has emerged as a powerful tool for studying hadronic charm 
decay, providing a new handle in dealing with FSI. In this context the recent high- 
statistics data have already furnished new interesting results: fully coherent analyses 
are now available for the KKn, tttttt and Ktttt final states of the D mesons. These 
decays are found to be dominated by two-body resonances with the unique excep- 
tion of D + — > K~7t + tt + . Comprehension of the whole decay pattern depends on 
a correct parametrization of the non-resonant contribution, which could influence 
evaluation of amplitudes and phases. Attempts to improve the non-resonant flat 
dynamics, allowing a variation over the Dalitz plot, have been proposed 0) and 
will be implemented in the next analyses. Looking to the future prospects: exciting 
data are expected from the new E831 run (now in progress) while E791 should give 
results pretty soon. When the quality of the data allows a study of CP violation in 
the charm-meson sector through the measurement of the D and D phases, a reliable 
estimate will require the ability to factor out the FSI induced phases. 
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Figure 7: Invariant tt + tt~ mass projections for the D + : a) low, b) high, c) sum; and 
.D+ : d) low, e) high, f) sum. The crosses are data, the continuous line is the fit 
result and the shaded area represents the background. 



